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ABSTRACT 

Context. Observations of protostellar disks indicate the presence of the magnetic field of thermal (or superthermal) strength. In such 
a strong magnetic field, many MHD instabilities responsible for turbulent transport of the angular momentum are suppressed. 
Aims. We consider the shear-driven instability that can occur in protostellar disks even if the field is superthermal. 
Methods. This instability is caused by the combined influence of shear and compressibility in a magnetized gas and can be an efficient 
mechanism to generate turbulence in disks. 

Results. The typical growth time is of the order of several rotation periods. 
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1. Introduction 

Protostellar disks require sufficiently strong turbulence to en- 
hance the efficiency of angular momentum transport. The origin 
of turbulence is often attributed to hydrodynamic and hydromag- 
netic instabilities that can arise in differentially-rotating, strati- 
fied gaseous disks. One of the candidates is magnetorotational 
instability (MRI), which can operate in a conductive flow if the 
angular velocity decreases with the cylindrical radius and the 
magnetic field is not strong (Velikhov 1959). The MRI has been 
studied in depth in the case of stellar and accretion disk con- 
ditions (see, e.g., Fricke 1969; Safronov 1969; Acheson 1978, 
1979; Balbus & Hawley 1991; Kaisig et al. 1992; Kumar et 
al. 1994; Zhang et al. 1994). Simulations of the MRI in disks 
(Hawley et al. 1995, Matsumoto & Tajima 1995, Brandenburg 
et al. 1995, Torkelsson et al. 1996, Arlt & Rudiger 2001) show 
the turbulence generated can significantly enhance the angular 
momentum transport. 

Most likely, however, the number of instabilities that can 
arise in astrophysical disks is quite large. An analysis of MHD 
modes in stratified accretion disks demonstrates a wide vari- 
ety of instabilities even in the case of simple magnetic geom- 
etry (Keppens, Casse & Goedbloed 2002). Therefore, the cur- 
rent point of view on the origin of turbulence in disks is likely 
highly simplified. Even pure, hydrodynamic origin of turbu- 
lence cannot be excluded (see Urpin 2003, Arlt & Urpin 2004, 
Dubrulle et al. 2005) despite the most efficient local linear hy- 
drodynamic instabilities advocated to date are not sufficiently 
efficient. Global instabilities, such as the baroclinic-like instabil- 
ity of Klahr & Bodenheimer (2003), are sensitive to boundary 
conditions (Johnson and Gammie 2006) and, therefore, are un- 
likely to drive turbulence in disks. Note, however, that such fac- 
tors as fast radiative cooling, high thermal diffusion, and large 
radial temperature gradients can strengthen the baroclinic feed- 
back and make the baroclinic instability more efficient (Petersen 



et al. 2007). Astrophysical disks are stratified and stratification 
can change stability properties of shear flows, providing either a 
stabilizing or destabilizing effect, depending on details of the 
disk structure. Recently, Dubrulle et al. (2005) have consid- 
ered one of the possible "strato-rotational" instabilities arising 
in the presence of both differential rotation and stable stratifi- 
cation. However, Brandenburg & Rudiger (2005) demonstrated 
that the growth rate of this instability decreases with an increas- 
ing Reynolds number, rendering the instability less relevant for 
astrophysical applications. Note also that convection can drive 
turbulence in disks with unstable stratification, but it induces in- 
ward turbulent transport of the angular momentum instead of the 
required outward one (see, e.g.. Stone & Balbus 1996). 

As it has been argued by Lesur & Longaretti (2005), nonlin- 
ear hydrodynamic instability that often occurs in linearly-stable 
flows at sufficiently large Reynolds numbers is likely inefficient 
in disks even for the most optimistic extrapolations of numerical 
data. At least, the subcritical transition to hydrodynamic turbu- 
lence cannot occur in quasi-keplerian flows at Reynolds num- 
bers up to ~ 10* (Ji et al. 2006). Of course, Reynolds numbers 
are much higher in real disks; nevertheless it seems that there is 
little room for nonlinear instability. 

The stability properties of protostellar disks are much differ- 
ent from those of accretion disks. The magnetic Reynolds num- 
ber is likely not very large in cold- and dense-protostellar disks 
because of a low electrical conductivity, and so the field cannot 
be treated as "frozen" into the gas (Gammie 1996). The effect 
of Ohmic dissipation on the MRI has been considered in the lin- 
ear (Jin 1996) and nonlinear regimes (Sano, Inutsuka & Miyama 
1998; Drecker, Rudiger & HoUerbach 2000). Fleming & Stone 
(2003) and Turner et al. (2006) treated turbulent mixing caused 
by the MRI in protostellar disks, taking into account magnetic 
diffusivity. They found that the midplane is shielded from cos- 
mic rays, and that MRI does not occur even under the most fa- 
vorable conditions. Nevertheless, turbulence can mix a fraction 
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of the weakly-conducting surface material into the interior, pro- 
viding some coupling of the midplane gas to the magnetic field. 
As a result, weak stresses can appear in the disk midplane. 

As first pointed out by Wardle (1999), poorly conducting 
protostellar disks can be strongly magnetized if electrons are the 
main charge carriers. Therefore, transport must be anisotropic 
with substantially different properties along and across the mag- 
netic field. In strongly magnetized gas the Hall component pro- 
vides the main contribution to the resistivity tensor, which pro- 
duces the electric field perpendicular to both the magnetic and 
electric current. The stability analysis by Wardle (1999) shows 
that the Hall eff'ect can provide either stabilizing or destabiliz- 
ing influences depending on the direction of the field. Balbus & 
Terquem (2001) have conducted a more general study of the role 
of the Hall term on MRI. They found that the HaU efifect quali- 
tatively changes the stability properties of protostellar disks and 
can lead to instability even if the angular velocity increases out- 
ward. These authors, however, did not take into account the ef- 
fect of gravity that is crucial for disks. Urpin & Rtidiger (2005) 
considered MRI under the combined influence of the Hall efl'ect 
and gravity, and also derived also the criteria of several other 
instabilities that can occur in protostellar disks. Sano and Stone 
(2002) investigated the effect of the Hall term on the evolution 
of the MRI in weakly-ionized disks using local axisymmetric 
simulations. These authors concluded the properties of the MRI 
depend essentially on the direction of the magnetic field as it 
is anticipated from the dispersion equation in a linear stability 
analysis. Salmeron & Wardle (2005) have also considered the 
properties of the MRI modified by the Hall effect. These authors 
argued that the MRI is active in protoplanetary disks over a wide 
range of field strengths and fluid conditions. The Hall conductiv- 
ity results in a faster growth of perturbations and extends the re- 
gion of instability. Recently, Livertz et al. (2007) and Shtemler et 
al. (2007) have considered the Hall MRI in a non-axisymmetric 
case. This type of instability is proposed as a viable mechanism 
for the azimuthal fragmentation of the protoplanetary disks and 
planet formation. The non-axisymmetric instability is caused by 
the combined effect of the radial stratification and HaU electric 
field. Note that the MRI and its modifications can be completely 
suppressed if the magnetic field is sufficiently strong (see Urpin 
1996, Kitchatinov & Rudiger 1997). 

If rotation is cylindrical and Q = f2(j), where Q. is the angu- 
lar velocity and s is the cylindrical radius, the critical magnetic 
field that suppresses the MRI is given approximately by the con- 
dition Ca/H ~ sQ' ~ Q where ca - B/ -sjA-np is the Alfven ve- 
locity. Since Q ~ c JH in the standard disk, where is the sound 
speed, we find that the MRI is suppressed if the magnetic field 
is superthermal, ca > Cj. Recent measurements (Hutawarakom 
& Cohen 1999, 2005, Donati et al. 2005) indicate that the mag- 
netic field can be strong in protostellar disks. This particularly 
concerns the innermost regions where the field strength reaches 
~ 1 kG at the radius ~ 0.05 AU (Donati et al. 2005). The field 
configuration includes the azimuthal component whose direction 
agrees with the radial field sheared by the disk differential rota- 
tion. The derived ratio of the azimuthal and radial fields in the 
surface layers is ~ 0.5. Note that the azimuthal component can 
be substantially stronger than the radial one in the disk interi- 
ors because dynamo theories (such as the o- - Q dynamo) usu- 
ally predict that the generated field should have a strong toroidal 
component. According to estimates by Donati et al. (2005), the 
equipartition field strength with roughly equal thermal and mag- 
netic pressures should be ~ 10^ Gauss at the radius ~ 0.05 AU 
in protostellar disks. This indicates the detected field is approx- 



imately of the thermal or even superthermal strength and can 
substantially modify or even suppress the MRI. 

Recently, Bonanno & Urpin (2006, 2007) have argued that a 
compressible, differentially-rotating flow is unstable if the mag- 
netic field has a non- vanishing, radial component. A remarkable 
feature of this shear-driven instability is that it can arise even 
if the magnetic field is superthermal. Bonanno & Urpin (2006, 
2007) have considered the instability in the case of a highly 
conductive plasma when the magnetic diffusivity plays no role. 
However, the instability can also arise in weakly-ionized pro- 
tostellar disks where Ihe magnetic field has a radial component 
detected in observations (Donati et al. (2005). The effect of a fi- 
nite diffusivity can be important for the onset of instability and 
for the properties of generated turbulence, particularly near the 
the mid-plane where conductivity is extremely low. In this pa- 
per, we consider the shear-driven instability in the conditions 
of protostellar disks and show that the conditions of instabil- 
ity in weakly and highly ionized plasmas differ substantially. 
Nevertheless, the shear-driven instability can manifest itself even 
in weakly-ionized protostellar disks. 

The paper is organized as follows. In Section 2, we con- 
sider the basic equations governing instabiUty in compressible 
dissipative fluids and we then derive the dispersion relation. In 
Section 3, we derive the criteria of instability and discuss condi- 
tions under which instability can occur in protostellar disks. The 
growth rate of instability is calculated in Section 4. Finally, in 
Section 5, we discuss the results obtained. 

2. Basic equations and dispersion relation 

The electrical conductivity is low in protostellar disks- and the 
magnetic field cannot be considered as "frozen" into the gas. 
The magnetic diffusivity is given by 77 = c^me/4ne^neT where 
me and are the mass and number density of electrons, respec- 
tively, and T is their relaxation time (see, e.g., Spitzer 1978). In 
protostellar disks, t is determined by the scattering of electrons 
on neutrals, then t = 1 /n(crv) where (crv) is the average product 
of the cross-section and velocity, and n is the number density of 
neutrals. Using the fitting expression for (crv) obtained by Draine 
et al. (1983), we have 

T] = 234 X 10^ f-^T^'^ cm^s-\ (1) 

where / = ng/n is the ionization fraction, and T2 = T/lOO K, 
with T being the temperature. 

The magnetic diffusivity can be anisotropic in some regions 
of protostellar disks (Wardle 1999) because the electron gas is 
magnetized despite a low temperature. The effect of the mag- 
netic field on transport properties is usually characterized by the 
magnetization parameter a^iB) = ojbT, where (Ob = eB/nieC is 
the gyrofrequency of electrons (see, e.g., Spitzer 1978). Again, 
using the fitting formula by Draine et al. (1983), we can estimate 
the magnetization parameter as 

a,(B) * 21 BuJIt;"^, (2) 

where the magnetic field B is measured in Gauss and nu = 
n/lO^'^cm-^. If fle(B) > 1, i.e. 

B > 0.048 ni4 ^/^\ Gauss, (3) 

then the electron transport is anisotropic and the magnetic diffu- 
sivity is represented by a tensor. In plasma of protostellar disks, 
the difference between components of the magnetic diffusivity 
parallel and perpendicular to the magnetic field is small (see. 
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e.g., Balbus & Terquem 2001). The Hall component is described 

as 



ae(B)T] = 



cB 

Anerii. 



(4) 



We work in cylindrical coordinates {s, ip, z) with the unit 
vectors (e,, e^, e^). The equations of compressible MHD read: 

v + (vV)v = -— +g+--^(VxB)xB, (5) 
p 4np 



p + V ■ (pv) = 0, 

p + v ■Vp + ypV • V = 0, 



B - V X (v X B) + V X 

4ne 



— (V xB)xB 

Me 



+7]V X (V X B) = 0, 



V • B = 0. 



(6) 
(7) 



(8) 
(9) 



Our notation is as follows: p and v are gas density and veloc- 
ity, respectively; p is pressure; y is the adiabatic index; and g is 
gravity. The third term on the l.h.s. of Eq. (8) describes the Hall 
effect. We cannot neglect this term because the magnetization 
parameter is large in protostellar disks. However, in what fol- 
lows, we will consider a special type of perturbation for which 
the Hall effect is unimportant. 

The basic state on which the stability analysis is performed 
is assumed to be quasi-stationary with the angular velocity Q - 
£lis) and B i= 0. We assume hydrostatic equihbrium for basic 
state: 



Vp 1 

— = D+ (V X B) X B 

p 4np 



D = g + Q^s. 



(10) 



We consider magnetic configurations where both the radial and 
azimuthal field components are present. The presence of a radial 
magnetic field and differential rotation leads to the development 
of the azimuthal field. In protoplanetary disks, however, the rate 
of stretching of the azimuthal field from is reduced by ohmic 
dissipation and Hall effect. Then, Eq. (8) yields the following 
steady-state condition: 



[Vx(vxB)]^ =.^Vx 



-(yxB)xi 



"I 



4nene 

+77[Vx(VxB)V (11) 
Assuming that > Bs, Eq. (11) can be transformed into: 



c / d iB^dBA d 
4ne \ds \ ne dz ) dz 



B^ IdB. 



ds 



Ba 



1 



77 A - — = sQ!Bs 



(12) 



If the Hall effect is negligible, the first term on the l.h.s. of this 
equation is small, and the generated toroidal field is stronger than 
the radial field by a factor of the order of the magnetic Reynolds 

number. Re,,,, 



B,D 



Re^Bv 



(13) 



where Re„ = s^I'H^/tj ~ CIH^/t] with H being the half- 
thickness of the disk (we assume sQ,' ~ Q). If the magnetic field 



is strong and the Hall effect dominates ohmic dissipation, then 
one estimates from Eq. (12) that. 



Bio 



Re 



1/2 



aJHB,) 



(14) 



In strongly magnetized disks, the generated toroidal field can 
be substantially weaker than what follows from the simplest es- 
timate (13). Generally, even if Re^ » 1, the toroidal field is 
comparable to the radial one if ~ Re„ despite a strong differ- 
ential rotation. It is possible, for example, that the toroidal and 
poloidal fields measured in FU Orionis are comparable because 
of the Hall effect. 

We consider the stability of axisymmetric short-wavelength 
perturbations. Small perturbations will be indicated by subscript 
1, while unperturbed quantities will have no subscript. The lin- 
earized momentum equation is: 



vi + (vi • V)v -I- (V • V)vi = -— -I- — 

P P 



p 



—(V xB)xB 
4np 



1 



-I- (VxBi)xB 

4np 

+—iVxB)xBi. 
4np 



(15) 



Taking into account that the unperturbed motion is rotation (v = 
sQe^) and using Eq. (10), we can transform this equation into 



vi -I- 2Q X vi -I- e^sQ'vis = — — + —D 

P P 

1 

+ [(V X Bi) X B -H (V X B) X Bi]. 

4np 



(16) 



The term proportional to pi/p is typically small in a short- 
wavelength approximation. Indeed, we can estimate pi ~ c^pi 
where c, is the sound speed. Then, the pressure term in Eq. (16) 
is (cj//l)(pi/p) where A is the lengthscale of perturbations. 
Observations indicate that the angular velocity can be smaller 
than the Keplerian angular velocity at least in a fraction of the 
disk volume (Donati et al. 2005) and, hence, D ~ g. Note that 
even departures by a factor of 2-3 from the Keplerian velocity re- 
ported by Donati et al. (2005) can produce a substantial increase 
in the local disk thickness compared to the "canonical" value 
~ QAs. Comparing the first two terms on the r.h.s. of Eq. (16), 
we conclude that the pressure term is larger if Cj//i » D ~ or 
A <K Cj/g. Since c]lg is the pressure scaleheight, this inequality 
is equivalent to the condition of applicabiUty of a short wave- 
length approximation. Therefore, we can neglect the term pro- 
portional to D in what follows. An additional simplification is 
related to the last term on the r.h.s. of Eq. (16) (the "curvature 
term" according to Pessah & Psaltis 2005). The toroidal field can 
often be stronger than the poloidal one and, as a result, the influ- 
ence of unperturbed electric currents can be important. However, 
if the toroidal field satisfies the condition: 



< {s/A)maxiBs,B^), 



(17) 



the contribution of "the curvature term" is small and can be ne- 
glected. 

Following similar transformations with the remaining 
Eqs. (6)-(9), we arrived at the linearized equations needed for 
stabiUty analysis. We consider perturbations with the space-time 
dependence oc exp(cr/ - ik ■ r), where k = {ks,0,k^ is the 
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wavevector. Then, the linearized MHD-equations read with ac- 
curacy in the lowest order in Als: 



oTi -I- 2n X Vi -I- BtpSDlvu 



ikpi 



-—{kxBOxB, 

api - ip{k ■ vi) = 0, 
crpi - iyp{k ■ Vi) = 0, 

(cr + (i),j)Bi = etpsQ.'Bis - i(B ■ k)vi + iB(k ■ vi) 

cik ■ B) , 



Anerit 



-kxB, 



k-B,^0, 



(18) 

(19) 
(20) 

(21) 

(22) 



where co,, = rjk^. 

The dispersion equation, corresponding to Eqs. (18)-(22) is 
rather complex in the general case. Therefore, we consider only 
a particular case when the wavevector is perpendicular to the 
magnetic field, k - B - 0. After some algebra, Eqs. (18)-(22) can 
be combined into a sixth-order dispersion relation, 

cr^ + a^cr' + a4<T^ + 030"^ + a2<J^ + a\cr + oq — 0. (23) 

The coefficients of this equation are expressed in terms of char- 
acteristic frequencies, 

as = 2a),j, 04 = 0)^ + col + 
a2 = +K^)+ iJ.K^(ol, 

where /! = l^/k^. The characteristic frequences are given by 

(Ajj„ = c\k^ , Ca= B/ ^J47Tp , = k^CA^CA,s^l' , 



where = ^fypJp is the sound speed. Eq. (23) describes six 
modes that can generally exist in a compressible, rotating, mag- 
netized gas. 

In the non-dissipative, incompressible, limit when the sound 
speed C5 is very large, Eq. (23) reduces to the dispersion relation 
for the inertial waves. 



0-2(0-2 -I- HK^) = 0. 



(24) 



This equation allows for unstable solutions if the Rayleigh crite- 
rion is fulfilled, <0. Note that Eq. (23) does not describe the 
magnetorotational instability since it does not occur for pertur- 
bations with k ■ B = 0. 

In non-dissipative linut, // — > 0, Eq. (23) recovers the disper- 
sion relation derived by Bonanno & Urpin (2006) for perturba- 
tions with k ■ B = 0, 



c 322 23 22 23 

cr + cr^ (lOq + k ) + cr oJg^-^ + crjiK cOq + /J.K <jJ^q = 0. 



(25) 



The authors argued that Eq. (25) has unstable solutions even 
when the criteria of the magnetorotational and Rayleigh instabil- 
ity are not satisfied. The condition of instabiUty is (Obq + and, 
hence, the shear-driven instability arises if 0. The paper by 



Bonanno & Urpin (2007) considers the dispersion relation for 
any wavevectors, but does so in a high-conductivity limit // — > 0. 

Eq. (23) can be simplified in many cases of interest if we 
take into account that, likely, the i-component of the magnetic 
field in disks is greater than the z-component, Bj » B^. For 
perturbations with ^ • B = 0, we have ks = -k^BJBg. Then, 



Bi 



Bl + Bl 



1. 



(26) 



Substituting ^ ^ \ into Eq. (23), we can transform it into 

(o-2 -1- i^)[cr'^ + Icjrjcr' + (w^ + u)l)cr^ + 

{oj^oj] + oj^ojI + o}l^)a- + ojIoj]-] = 0. (27) 

Two roots of this equation describe the inertial waves which can 
be unstable only if the Rayleigh criterion is satisfied. Other four 
modes are described by the dispersion relation 

-I- 2(0r,cr^ + {a^ + o^)(t^ + (iOrjoJl + oJrjO^ + co\^cr 

+0)^0)1 = 0. (28) 

This equation describes fast and slow magnetoacoustic waves, 
and we consider the stability of these modes. Note that simpli- 
fication (26) is made for mathematical convenience rather than 
for physical relevance. In fact, the dispersion relation (28) for 
fast and slow magnetosonic waves will change little if /i 9^ 1. 
This particularly concerns the modes with |o-| > Q because the 
dispersion equation for them will differ from Eq. (28) only by 
terms of the order of (fi - \)i^lo^ in coefiicients. Therefore, our 
results can be applied with a sufficient accuracy also for the per- 
turbations with /I 5t 1 . 



3. Criteria of instability in protostellar disks 

The Hurwitz theorem states that an equation of the fourth order, 

a'^ + bscr^ + bicr^ + hcr^ +bo = 0, (29) 

has at least one root with a positive real part (unstable mode) if 
one of the following inequalities 



Z73 < , bQ<0, 
b-ib2 - fei < , 
bi{bib2-bi)-b\bo<0 , 



(30) 
(31) 
(32) 



is fulfilled (see Aleksandrov et al. 1985). Conditions (30) are 
never satisfied in protostellar disks. Therefore, the instability 
arises if one of the inequahties (31) or (32) are satisfied. These 
inequalities can be rewritten for Eq. (28) as 



,2, ,2 , 



<0, 



(33) 
(34) 



Generally, the expressions on the l.h.s. of both these inequalities 
can have a positive or negative sign, depending on the value of 
Wfin- Therefore, Eq. (33) and (34) impose restrictions on the rate 
of differential rotation which should be greater than some critical 
value. Note that, if the azimuthal field is generated by winding 
up the radial field, then BjB^Q' > and, hence, o?g^ > as well. 
Therefore, we consider only this case. Denoting 6 = CasCaipIc\, 
we have from Eq. (33) 



HsD! 

6 > 

ca caH 



(kHT + {kHf 



(35) 
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Fig. 1. The region of parameters where the instabihty can occur 
in accordance with Eq. (32) for different values 6. 



Fig. 2. The region of parameters where instability can occur in 
accordance with Eq. (33). The values of {6,P\ are [-0.5, 1] (solid 
lines), [-0.1, 1] (dashed Unes), and [-0.1, 10] (dashed-and-dotted 
lines). 



This condition can be satisfied only if the differential rotation 
is sufficiently strong. Since Q' < in astrophysical disks, we 
consider only this case. Then, we have 5 <0 because BsB^Qf oc 
6Q.' > 0. The r.h.s. of Eq. (35) is a monotonously increasing 
function of kH and is a minimum for the minimal possible value 
of kH. Since a short-wavelength approximation applies only if 
kH > 1, we can estimate the domain of instability supposing 
kH = 1 in inequality (35). Then, in order to fulfill the instability 
condition (35) for short-wavelength perturbations with kH > 1, 
the rotation shear has to satisfy, at least, the inequahty: 



HsQ.' T] 

6 > — 

ca caH 



+ 1 



(36) 



In Fig. 1, we plot the regions of instability given by Eq. (36) for 
different negative values of 6. At given 6, the instability can oc- 
cur if the value of \Hs^' /ca\ is above the corresponding curve 
(note that Q' is negative in Fig. 1). If the magnetic and thermal 
energies are comparable and departures from the Keplerian disk 
are not large, then ca/H ~ c^/H ~ £1, and HsQ.'/Cs ~ sQ.'/Q. « 
3/2. Inequality (36) can be satisfied only for relatively small val- 
ues of TjjHcA < or \6\ ~ Bg/B^ > l/Re™. This condition can 
be satisfied only if the Hall parameter is large > I, because 
6 ~ l/Rsm (see Eq. (13)) in the case < 1. Condition (36) can 
also be fulfilled in a non-Keplerian disk. For a non-Keplerian 
disk, the quantity HsQ.' /ca can vary in a wider range even if 
the magnetic and thermal energies are comparable, but condi- 
tion (36) still cannot be satisfied for relatively large values of 
r]/HcA- 

By analogy, we can also transform condition (34). The min- 
imal unstable rotation shear corresponds again to the maximal 
possible wavelength. Again, the maximum wavelength can be 
estimated from the condition kH ~ I. For such wavevectors. 



Eq. (36) can be transformed into: 



,(HsD. 



\ Ca 



I \ CA }Hca\h^cI ^) 



H^cir 



(37) 



where jS = c^/c^. 

In Fig. 2, we show the region of parameters where the in- 
stability can occur in accordance with Eq. (37). Inequality (37) 
is quadratic in HsQ.' /ca and can be satisfied if this parameter is 
larger/smaller than the largest/smallest root of the corresponding 
quadratic equation. Therefore, for a given combination of yS and 
S, the regions of instabihty lie above the upper line and below 
the lower line. 



4. Growth rate of instability 

To calculate the instability growth rate, it is convenient to intro- 
duce dimensionless quantities. 



n ^ _Cs_ 
DB^ ' ^ Hn' ^ 

. Ca^Cas 



= kH , 



(we assume BsB^Q! > 0). Then, Eq. (24) transforms into 

+ lax^T^ + +T^x\(?}? +e\\ + IIP)] + 
+reV[ax\2 + qdip] + a^s^x^ = 0. 



(38) 



The dependence on the wavelength is characterized by the pa- 
rameter in this equation. The parameter a is of the order of 
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Fig. 3. The dependence of the growth rate on x = kH for the Fig. 4. The same as in Fig. 3 but for a nonKeplerian disk with 
Keplerian disk, with q = 3/2 and e = 1, and for yS = 0.1 (soUd q = I and e = 3. 
line), 1 (dashed line), and 10 (dashed-and-dotted Une). The other 
parameters are J = 0.1 and a = 0.001. 



0.08 



1 /Re^. We solved eq. (38) numerically for different values of the 
parameters by computing the eigenvalues of the matrix whose 
characteristic polynomial is given by Eq. (38) (see Press et al. 
1992 for details). 

In Fig. 3, we plot the dependence of the growth rate on x 
for the Keplerian disk with q - 312 and e = 1 and for sev- 
eral values of the parameter p. In all cases, the growth rate de- 
creases monotonically with decreasing wavelength A = 2njk be- 
cause ohmic dissipation is more efficient for perturbations with 
a shorter wavelength. For any ratio of the magnetic and gas pres- 
sures, the instability does not occur if kH > 12. If the magnetic 
pressure is greater than the gas pressure (J3 < I), the growth rate 
can reach relatively large values F 0.04 - 0.06. This corre- 
sponds to the growth time only 3-4 times longer than the ro- 
tation period. The magnetorotational instability cannot arise in 
such a strong magnetic field no matter how large is the electrical 
conductivity of gas. Even if k is not perpendicular to B, the MRI 
can occur only if the magnetic field and wavevector satisfy the 
condition k~c^ < 2sQ.Sl'. Assuming that sCl' ~ Q, we can trans- 
form this inequality approximately into Q > c^k ~ xcJHy/p. 
Since Cs/H ~ Q. in the Keplerian disk, we find that MRI arises if 
xl ^Ip < 1 . In the case yS < 1 , this condition never applies and the 
MRI is suppressed. On the contrary, the shear-driven instability 
turns out to be rather efficient in such a strong field. 

In Fig. 4, we plot the same dependence, but for ffie nonKe- 
plerian disk wiffi ^ = 1 and e = 3. The offier parameters are 
same. Qualitatively, the dependences are very similar to ffiose 
shown in Fig. 3. The instability grows faster in the disk with 
larger /?, which seems to be a general feature of the shear-driven 
instability. However, ffie growth rate is typically smaller for a 
non- Keplerian disk. The range of unstable wavevectors is lim- 
ited approximately by kH < 10. Note that ffie half-thickness H 
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Fig. 5. The dependence of the growffi rate on yS for ffie Keplerian 
disk wiffi q = 3/2, e=\,6 = 0.1, and a = 0.001. The soUd and 
dashed lines correspond to x = 5 and 10, respectively. 



is essentially larger in a non-Keplerian disk than in the Keplerian 
one. 

The dependence of ffie growth rate on p'^ = c\lcl is shown 
in Fig. 5 for the disk with q - 3/2 and s - 1. At a given wave- 
length, the growth rate increases with 1 //? slowly approaching a 
saturation value that depends on q. The saturation growth rate is 
~ 0.06Q and 0.025Q for ^ = 5 and 10, respectively. This corre- 
spond to the growffi time of about ~ 3 and ~ 6 rotation periods. 
The dependence of the shear-driven instability on the magnetic 
field is unusual, at least, in ffie domain of parameters because 
usually a strong magnetic field tends to suppress other types of 
instabiUty. Fig. 6 shows ffie dependence of ffie growth rate on the 
parameter a, which is approximately equal to the inverse mag- 
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Fig. 6. The dependence of the growth rate on a for q - 312, e - 
1,5 -Q. I, and yS = 0.1. The soHd and dashed curves correspond 
to ^ = 5 and 10, respectively. 



netic Reynolds number, a ~ 1/Re„,. Therefore, higher values 
of a correspond to a higher dissipation rate, and the instability 
should be suppressed at large a. Indeed, the growth rate goes to 
if a ~ 0.001 - 0.01, depending on the wavelength. Hence, one 
can expect that the shear-driven instability occurs even in disks 
with a relatively small magnetic Reynolds number ~ 100. 

5. Discussion 

We have considered the new instability that can arise in proto- 
stellar disks under the combined influence of the magnetic field, 
differential rotation, and compressibility. To illustrate the main 
qualitative features of the instability, we analyzed a particular 
case of perturbations with the wavevector k perpendicular to the 
magnetic field B. In this case, the standard MRI does not occur 
because its growth rate is oc ■ B). The considered instability is 
related to shear and compressibility of a magnetized gas. In the 
incompressible limit (c, — > oo) we have, from Eq. (27), 

(cr2 + K^)(o- + tj^f = 0, (39) 

and the instability does not occur. This is a major difference 
from other instabilities caused by differential rotation, such as 
the Rayleigh or magnetorotational instabilities that can occur in 
incompressible limit. 

The magnetic field is likely thermal or superthermal in pro- 
tostellar disks (see Donati et al. 2005). The considered insta- 
bility can arise even in a very strong magnetic field when the 
MRI is suppressed completely. It is known that the MRI does 
not occur if the magnetic field satisfies the condition B > B^r - 
(sD.' /k) ^J4■np. The shear-driven instability can operate even if 
the field is stronger than B^r- The growth time of instability is 
still rather short (several rotation periods) in the gas where the 
magnetic pressure is greater than the thermal pressure. 



The criteria for the considered instability can be satisfied 
in protostellar disks. Likely, all disks have both radial and az- 
imuthal magnetic fields. These magnetic fields are found in ob- 
servations, as well as in numerical simulations. Superthermal 
magnetic fields do not suppress the instability crucial in cold pro- 
tostellar disks. The instability can arise even in a low conductive 
gas where the magnetic Reynolds number associated with rota- 
tion is relatively small, Re,„ ~ 100. 

Our paper considers only a local axisymmetric instability. 
Most likely, the same mechanism of instability will efficiently 
destabilize perturbations with the wavelength comparable to the 
disk scale height. Dissipative effects are less important for these 
perturbations, and the unstable modes can grow even faster A 
global instability caused by a combined influence of differential 
rotation and compressibility will be considered elsewhere. The 
turbulence that could be generated by the considered instabil- 
ity can be strongly anisotropic in the (i,z)-plane because both 
the criteria and growth rate are sensitive to the direction of the 
wave vector. The generated turbulence, however, might be ef- 
ficient in the radial transport of angular momentum in strongly 
magnetized protostellar disks. 
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